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A sequence of moments obtained from statistical trials encodes a classical probability distribution. 
However, it is well-known that an incompatible set of moments arise in the quantum scenario, 
when correlation outcomes associated with measurements on spatially separated entangled states are 
considered. This feature viz., the incompatibility of moments with a joint probability distribution is 
reflected in the violation of Bell inequalities. Here, we focus on sequential measurements on a single 
quantum system and investigate if moments and joint probabilities are compatible with each other. 
By considering sequential measurement of a dichotomic dynamical observable at three different time 
intervals, we explicitly demonstrate that the moments and the probabilities are inconsistent with 
each other. Experimental results using a nuclear magnetic resonance (NMR) system are reported 
here to corroborate these theoretical observations viz., the incompatibility of the three-time joint 
probabilties with those extracted from the moment sequence when sequential measurements on a 
single qubit system are considered. 
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I. INTRODUCTION 



The issue of determining a probability distribution 
uniquely in terms of its moment sequence - known as 
classical moment problem, - has been developed for more 
than 100 years [1, 2]. In the case of discrete distributions 
with the associated random variables taking finite values, 
moments faithfully capture the essence of the probabili- 
ties i.e., the probability distribution is moment determi- 
nate [3]. 

In the special case of classical random vari- 
ables Xi assuming dichotomic values Xi — ±1, 
it is easy to see that the sequence of mo- 
ments [4] Mnin....n. = {X^'X1^-...XI 



Ii'nin2...nk,ni,n2, ... rife = 0, 1 as, 
P{xi,X2,...,Xk) ^ :^ ^ 

ni,...nk- 



2k Z^ •"! -^2 

ni,...nfc=0,l 



n\ n2 Tik 

X-\ Xg . . . Xj^ f^ni ...rik 



1 



E 

ni,...,nfc=0,l 



(1) 



:r) 



X-y^ X^ ■■■ X^.'^ P(Xi,X2, ■■.,Xk) 



where 



Xl,X2,... Xk=±l 

ni,n2, . . . ,?T.fe = 0,1, can be readily inverted to 
obtain the joint probabilities P{xi,X2, ■ ■ ■ ,Xk) 
uniquely. More explicitly, the joint probabilities 
P{xi, X2, . . . , Xk) are given in terms of the 2*^ moments 
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Does this feature prevail in the quantum scenario? 
This results in a negative answer as it is wellknown that 
the moments associated with measurement outcomes 
on spatially separated parties are not compatible with 
the joint probability distribution. This feature reflects 
itself in the violation of Bell inequalities. In this 
paper we investigate whether moment-indeterminacy 
persists when we focus on sequential measurements 
on a single quantum system. We show that the dis- 
crete joint probabilities originating in the sequential 
measurement of a single qubit dichotomic observable 
X{ti) — Xi at different time intervals are not consistent 
with the ones reconstructed from the moments. More 
explicitly, considering sequential measurements of Xi, 
X2, X3, we reconstruct the trivariate joint probabil- 
ties Pfj^(xi,X2,X3) based on the set of eight moments 



{(Xi), (X2), (X3), (XiX2),(X2X3), (X1X3), (X1X2X3)} 

and prove that they do not agree with the three-time 
joint probabihties (TTJP) Pj_{xi,X2,X3) evaluated di- 
rectly based on the correlation outcomes in the sequential 
measurement of all the three observables. Interestingly, 
the moments and TTJP can be independently extracted 
experimentally in NMR system - demonstrating the 
difference between moment inverted three time proba- 
bilities with the ones directly drawn from experiment, 
in agreement with theory. For obtaining TTJP directly 
we use the procedure of Ref. [5] and for extracting 
moments we extend the Moussa protocol [6] to a set of 
non-commutating observables. The specifics are given in 
the experimental section. 

Disagreement between moment inverted joint proba- 
bilities with the ones based on measurement outcomes 
in turn reflects the inherent inconsistency that the fam- 
ily of all marginal probabilities do not arise from the 
grand joint probabilities. The non-existence of a le- 
gitimate grand joint probability distribution, consistent 
with the set of all pairwise marginals is attributed to be 
the common origin of a wide range of no-go theorems 
on non-contextuality, locality and macrorealism in the 
foundations of quantum theory [7-14] . The absence of a 
valid grand joint probability distribution in the sequen- 
tial measurement on a single quantum system is brought 
out here in terms of its mismatch with moment sequence. 

We organize the paper as follows. In Sec. II we begin 
with a discussion on moment inversion to obtain joint 
probabilities of three classical random variables assum- 
ing dichotomic values ±1. We proceed in Sec. Ill to study 
the quantum scenario with the help of a specific example 
of sequential measurements of dichotomic observable at 
three different times on a spin-1/2 system. We show that 



the TTJP constructed from eight moments do not agree 
with those originated from the measurement outcomes. 
Sec. IV is devoted to report experimental results with 
NMR implementation on an ensemble of spin-1/2 nuclei, 
demonstrating that moment constructed TTJP do not 
agree with those directly extracted. Section V has con- 
cluding remarks. 



II. RECONSTRUCTION OF JOINT 

PROBABILITY OF CLASSICAL DICHOTOMIC 

RANDOM VARIABLES FROM MOMENTS 

Let X denote a dichotomic random variable with out- 
comes X = ±1. The moments associated with statis- 
tical outcomes involving the variable X are given by 
^„ = (X") = E.=±ia^"-P(a^). ri = 0,1,2,3,..., where 
< P{x = ±1) < l;Ex=±i^(^) = 1 are the corre- 
sponding probabilities. Given the moments fiQ and fj,i 
from a statistical trial, one can readily obtain the prob- 
ability mass function: 

i.e., moments determine the probabilities uniquely. 

In the case of two dichotomic random variables 
Xi, X2, the moments /i„i,„2 = (X"^ ^2'^) = 

y x^^ X2^ P{xi,X2), ni,n2 — 0,1... encode 

Xi— ±1,X2— ±1 

the bivariate probabilities P{xi,X2)- Explicitly, 



Xi,X2—:tl 

Min = ^ xiP{xi,X2) ^ ^ xi P(xi), 

xi,X2—dil a:i— ±1 

= P(i, 1) + P(i, -1) - P(-i, 1) - P(-i, -1) 

^J■01 = ^ X2Pixi,X2) =^ y^ X2P{X2) 

X-i^X2 — ±l 2:2 

= P(l, 1) - P(l, -1) + P(-l, 1) - P(-l, -1) 

Mil = Y. 2:1x2 P(a;i,X2)=P(l,l)-P(l,-l)-P(-l,l)+P(-l,-l) 

a:i,a:2— ±1 



(2) 



r 



Note that the moments /iio, /^oi involve the marginal 
probabilities P{xi) = Yl,x^=±\ P{x\,X2), P{x2) = 
X^x =±1 ^(^11^2) respectively and they could be evalu- 
ated based on statistical trials drawn independently from 
the two random variables X\ and X2. 



Given the moments /J.oo,MiO:Moijy^ii the reconstruc- 



tion of the probabilities P(a;i,a;2) is straightforward: 



L ^2 

2^2 A^ni "2 



P{X1,X2) = - Yl ^l'" 

ni,n2— 0,1 

= 7 E ^r^2M^r^2"^)- (3) 



ni,n2—0,l 



Further, a reconstruction of trivariate joint probabiH- 
ties P{xi,X2,X3) requires the following set of eight mo- 
ments: {^000 = 1,M100 = (-^l), MOIO = {X2),^J-010 = 

{X3), 1^110 = (-''^i -''^2)! Moil — (-'^2 -'^a), Mioi = 
{X1X3), fiiii — {X1X2X3)}. The probabilities are re- 
trieved faithfully in terms of the eight moments as, 



•^2 '''3 r'ni "2 "3 



P(a;i,X2,X3) = - ^ a;"i 

ni,n2,n3— 0,1 

= 1 E :^r4^4M^r^2"^x 



n3\ 
3 / 



ni,n2,n3— 0,1 



(4) 



It is implicit that the moments /iioOiMoiOi A'ooi are 
determined through independent statistical trials in- 
volving the random variables Xi , X2 , X^ separately; 
A'iiOjMoiijMioi are obtained based on the correla- 
tion outcomes of (Xi,X2), {X2,X3) and {XijX^) 
respectively. More specifically, in the classical 

probability setting there is a tacit underlying as- 
sumption that the set of all marginal probabil- 
ities P{xi),P{x2), P{x3),P{xi,X2),Pix2, X3), P(xi, Xa) 

are consistent with the trivariate joint probabilities 
P(a;i, 2:2, X3). This underpinning does not get imprinted 
automatically in the quantum scenario. Suppose the ob- 
servables Xi , X2 , ^3 are non-commuting and we consider 
their sequential measurement. The moments /iioo = 
(-^i),A*oio = {X2), Mooi = (-^3) may be evaluated 
from the measurement outcomes of dichotomic observ- 
ables Xi , X2 , X3 independently; the correlated statistical 
outcomes in the sequential measurements of (^1, X2), 
{X2; ^3) and {Xi, X3) allow one to extract the set 
of moments ^no = (-^1-^2), Mon = (^2-^3), Mioi = 
{XiX^); further the moment /im — (Xi X2 X3) is eval- 
uated based on the correlation outcomes when all the 
three observables are measured sequentially. The joint 
probabilities P^{xi,X2,X3) retrieved from the moments 
as given in (4) differ from the ones evaluated directly in 
terms of the correlation outcomes in the sequential mea- 
surement of all the three observables . We illustrate this 
inconsistency appearing in the quantum setting in the 
next section. 



III. QUANTUM THREE-TIME JOINT 
PROBABILITIES AND MOMENT INVERSION 



Let us consider a spin-1/2 system, dynamical evolution 
of which is governed by the Hamiltonian 



H = -huax 



(5) 



We choose z-component of spin as our dynamical observ- 
able: 



X, = X(U) = a,{U) 

— Uz cosLoti + Oy sinwii. 



(6) 



where lJ{ti) — e *°"^"*i/2 — [/.^ and consider sequen- 
tial measurements of the observable Xi at three different 
times ti = 0, ^2 = Ai, ^3 = 2 Ai: 

Xi — dz 

X2 ~ ijz{l^t) = Oz cos((x>At) + cTj, sin(a;At) 

Xz = cTz{2At) = <7z cos{2ujAt) + Oy sin(2wAt). (7) 

Note that these three operators are not commuting in 
general. 

The moments (Xi), (X2), {X^) are readily evaluated to 
be 

/^lOO = (-^1) = Tr[/5inCr^] = 0, 

Molo = (X2)=Tr[/}i„(7,(Ai)]=0, 
Mooi = (X3)=Tr[pi„fT3(2Ai)]=0. 

when the system density matrix is prepared initially in 
a maximally mixed state pin = 1/2. The probabilities 
of outcomes Xi = ±1 in the completely random initial 
state are given by P{xi — ±1) — Tr[pi„n2,.] = i, where 

n^i = |a^i)(2;i| is the projection operator corresponding 
to measurement of the observable Xi . 

The two-time joint probabilities arising in the sequen- 
tial measurements of the observables Xi^Xj^ j > i are 
evaluated as follows. The measurement of the observable 
Xi yielding the outcome Xi = ±1 projects the the density 

operator to p^^ — ^'. '"- ^? . Further, a sequential mea- 
surement of Xj leads to the two-time joint probabilities 
as, 



P{Xi,Xj) = P{Xi) P{Xj\Xi) 

= Tr[n,,pinnx,n,j 



{Xilpin \Xi) \{Xi\Xj) 



(8) 



We evaluate the two-time joint probabilities associated 
with the sequential measurements of (Xi,X2), {X2,X3), 



and {Xi,Xs) explicitly: 

P{xi,X2) = - [1 + xi X2 cos{ujAt)] 
P{x2, X3) = - [1 + X2 X3 cos{ujAt)] 
P{xi,X3) = - [1 + xi X3 cos(2wAi)]. 
We then obtain two-time correlation moments as, 

^110 = (i^li'2) = ^ XiX2P{xi,X2) 

= cos(wAt) 

flOll = {X2X3) = ^ X2X3P{X2,X3) 

= COS (w At) 
^llOl^ {X1X3) = ^ XiX3P{xi,X3) 

X-i,X3—dzl 

= cos{2ujAt). 



(9) 
(10) 

(11) 



(12) 
(13) 
(14) 



Further, the three-time joint probabilities P{xi,X2,X3) 
arising in the sequential measurements of Xi,X2, fol- 
lowed by X3 are given by 



where p,, = "?f/^^. We obtain, 

p{xi,x2,x3) = Trpa;^ n^i piniia^j n^^^n^a] 

= {Xl\pin\xi) \{xi\x2)\'^ \{X2\X3)\'^ 
_ P{XI,X2)P{X2,X3) 

{X2\pin\x2) 
_ P{XI,X2)P{X2,X3) 
P{X2) 

where in the third line of (16) we have used (8). 



(16) 



The three-time correlation moment is evaluated to be, 

^111 = (i'lX2^3) = X! XiX2X3P{xi,X2,X3) 

2;i,2;2,K3=±l 

= 0. (17) 



P{xi,X2,X3) = P{xi) P{x2\xi) Pix3\xi,X2) 

= Tr[/5i„ n^J Tr[p^, U^^] Tr[px2 ^x^] 



(15) 



From the set of eight moments (8), (12) and (17), we 
construct the TTJP (see (4)) as. 



P^(l,l,l) = - [1 + 2 cos(wAt)-fcos(2LjAt)]=F^(-l, -1,-1), 
o 

^•(-1,1,1) = ^ [l-cos(2c^At)]= P^(-l, -1,1) = F^(l, 1,-1) = P^(l, -1,-1), 
P^(l,-l,l) = - [1-2 cos(a;At)+ cos(2wAt)] =P^(-1, 1,-1). 



(18) 



On the other hand, the three dichotomic variable quantum probabilities P{xi, X2, X3) evaluated directly are given by. 



Pd(l,l,l) = -[l + cos{uAt)f^Pd{-l,-l,-l), 
^d(-l,l,l) = ^[l-cos2(t^At)]=P<i(-l,-l,l) = Pd(l,l,-l) = Frf(l,-l,-l), 
Pd(l,-l,l) = i[l-cos(L.Ai)]' = P<i(-l,l,-l). 



(19) 



r 



Clearly, there is no agreement between the moment 
inverted TTJP (18) and the ones of (19) directly eval- 
uated. In other words, the TTJP realized in a se- 
quential measurement are not invertible in terms of the 
moments - which in turn reflects the incompatibility 
of the set of all marginal probabilities with the grand 
joint probabilities Pd{xi,X2,X3). In fact, it may be ex- 



plicitly verified that P(xi,a;3) ^ Z]i^2=±i Pdixi,X2,X3). 
Moment-indeterminacy points towards the absence of a 
valid grand probability distribution consistent with all 
the marginals. 

The TTJP and moments can be independently ex- 
tracted experimentally using NMR methods on an en- 
semble of spin-1/2 nuclei. The experimental approach 
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FIG. 1: Circuit for finding three-time probability. Grouped 
gates represent measurement in rotated basis and controlled 
gates can be CNOT or anti-CNOT as explained in the text. 



and results are reported in the next section. 



IV. EXPERIMENT 

The projection operators at time t = {Xi = Cz) 
are {II^^o = |2^i')(a;^|}a:0=o.i. This measurement basis is 
rotating under the unitary [7^, resulting in time depen- 
dent basis given by, n^.* = Ultl^oJJi. While doing ex- 
periments it is convenient to perform the measurement 
in the computational basis as compared to the time de- 
pendent basis. This can be done as follows: We can ex- 
pand the measurement on an instantaneous state p{ti) as, 



n^t/5(i,)n^ 



r7tfT 



[//n^o U.,p{U)Ul n^oC/,. Thus, mea- 



suring in time dependent basis is equivalent to evolving 
the state under the unitary Ui, followed by measuring 
in the computational basis and lastly evolving under the 
unitary Uj. 

The probabilities of measurement outcomes can be en- 
coded onto the ancilla qubits with the help of CNOT 
(or anti-CNOT) gate. To see this property consider a 
one qubit general state (for system) and an ancilla in the 
state |0)(0|, then the CNOT gate encodes the probabili- 
ties as follows 

(po|0)(0| +pi|l)(l| + a|l)(0| + at|0)(l|)^ ® |0)(0U 
4, CNOT 

|0)(0|s ®Po|0)(OU + |l)(l|s ®Pi|l)(lU 
+ |l)(0|s ^ a|l)(OU + |0)(l|s ^ a^|0)(lU. 

Now measuring the diagonal terms of the ancilla qubit, 
we can retrieve po and pi . 

We have employed a model as shown in Fig. 1 for 
measuring TTJP [5]. The grouped gates represent the 
measurements in the rotated bases. The controlled gates 
shown can be either CNOT or anti-CNOT gate. We re- 
quire both CNOT and anti-CNOT gate to perform the 
'ideal negative result measurement' (INRM) procedure to 
measure the TTJP noninvasively, as proposed by Knee 
et.al. [15]. The idea behind the INRM procedure is as 
follows: consider a gate which interacts with the ancilla 
qubit only when the system qubit is in state |1). By ap- 
plication of such a gate we can noninvasively obtain the 
probability of the measurement outcomes when the sys- 
tem was in |0) state. Similarly if we have a gate, which 



can interact with ancilla only if the system qubit is in |0) 
state then we can noninvasively obtain the probability of 
the measurement outcomes of the system state being in 
|1). These criteria are fulfilled by the CNOT gate and 
the anti-CNOT gate respectively. 

Circuit shown in Fig. 1 has two controlled gates for 
encoding the outcomes of first and second measurements 
on to the first and second ancilla qubits respectively. The 
third measurement need not be non-invasive since we are 
not concerned with the further time evolution of the sys- 
tem. A set of four experiments are to be performed, with 
following arrangement of first and second controlled gates 
for measurement of the TTJP: (i) CNOT; CNOT, (ii) 
anti-CNOT; CNOT, (iii) CNOT; anti-CNOT, and (iv) 
anti-CNOT; anti-CNOT. 

The propagators Ui = e~"^^"*'/^ is realized by the cas- 
cade HC/dH, where H is the Hadamard gate, and the 



delay propagator Ud = e 



-ia^ujti/2 



corresponds to the z- 



precession of the system qubit at w = 27rl00 rad/s reso- 
nance off-set. The diagonal tomography was performed 
at the end to determine the probabilities [5]. 

The three qubits were provided by the three ^^F nu- 
clear spins of trifluoroiodoethylene dissolved in acetone- 
D6. The structure of the molecule is shown in Fig. 2(a) 
and the chemical shifts and the scalar coupling values 
(in Hz) in Fig. 2(b). The effective ^^F spin-lattice (T^) 
and spin-spin (Ti) relaxation time constants were about 
0.8s and 6.3 s respectively. The experiments were car- 
ried out at an ambient temperature of 290 K on a 500 
MHz Bruker UltraShield NMR spectrometer. The first 
spin (Fi) is used as the system qubit and, other spins 
(F2 and F3) as the ancilla qubits. Initialization involved 
preparing the state, ^'^l + e{^\ls® |00)(00|a} where 
e ^ 10^^ is the purity factor [16]. The pulse sequence to 
prepare this state from the equilibrium state is shown in 
Fig. 2(c). All pulses were numerically optimized using 
the GRAPE technique [17] and had fidelities better than 
0.999. 

With our choice of measurement model (Fig. 1) we find 
a striking agreement with theoretical results on TTJP 
(19) . One might have also run the post measured state 
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FIG. 2: The molecular structure of trifluoroiodoethylene (a), 
and corresponding chemical shifts and J-coupling values(in 
Hz) (b), and the pulse sequence for the preparation of initial 
state (c). The open pulses are vr pulses and r = 1/(4J23). 



resulting after the first dashed block in Fig. 1 through an 
arbitrary CP map before the next step. However such 
post processing CP map would have affected the results. 
In other words, our measurement scheme provides an op- 
timal procedure to preserve the state information, thus 
resulting in an excellent agreement of experimental re- 
sults on TTJP with theoretical prediction (see Fig. 5). 

For calculating the moments we utilize the Moussa 
protocol [6], which requires only two spins in our 
case. We utilize Fi as the system and F2 as the 
ancilla qubit. F3 was decoupled using tt pulses and 
the initialization involved preparing the state, -i^l + 
e {ils (g) |+)(+U <^ |0)(0|}, which is obtained by apply- 
ing the Hadamard gate to F2 after the pulse sequence 
shown in Fig. 2(c). The circuit for measuring moments 
by Moussa protocol is shown in Fig. 3 and it proceeds as 
follows, 

p^l+X+l 

icXi 

pXl®|0)(l|-fXip®|l>(0| + 

p(E)\0){0\+XipX\(E)\l){l\ 

icX2 

pX\xl(g>\0){l\+X2X,p(g,\l){0\ + 

p (g) |0) (0| -I- X2XipX\xl (g) |1) (1| 

icXg 

pX|xJX^(8) |0)(1| +X3X2X1P0 |1)(0| + 



/5(g) 



+ X^X2XipX[XlXl(g,\l){l 



where, cXj represents the controlled gates and p is the 
initial state of the system. The state of the ancilla qubit 
(pa) at the end of the circuit is given by. 

Pa = \0){l\Tr{pXlxlxl) + \l){0\Tr{X3X2X,p) 

+ |0)(0|Tr(p) + \l){l\Tr{X3X2X,pXlxlxl). 

Moussa protocol was originally proposed for commutat- 
ing observables, however, it can be easily extended to 
non-commutating observables. The NMR measurements 
correspond to the expectation values of spin angular mo- 
mentum operators Ix or /j,[18]. The measurement of the 
Lj. for ancilla qubit at the end of the circuit gives: 

TT[pah] - Tr[X3l2^ip]/2 + TiipXlxlxl]/2. (20) 



X, 



X 



X 



i+x+i- 



X — X — lH— ^^ 



FIG. 3: Moussa Protocol for obtaining the 3-time correlated 
moments. One and two time moments can be calculated using 
the appropriate number of controlled gates. 




FIG. 4: Moments obtained experimentally from Moussa Pro- 
tocol. The symbols represent experimentally obtained values 
of the indicated moments with the solid lines showing the 
corresponding theoretical values. 
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FIG. 5: Three-time joint probabilities (TTJP): The solid 
curve represents the probabihties obtained directly and the 
dashed curve the probabihties obtained by inverting the mo- 
ments. The symbols represents the experimental data. 



If, Xi,X2,X3 commute, then the above expression gives 
Tr[/3XiX2X3]. In case of non-commuting hermitian ob- 
servables, we also measure expectation value of ly, which 
gives : 

iTvipaly] = Tt[X3X2Xip]/2 - TT[pXlxlxl]/2. (21) 

From (20) and (21) we can calculate Tt[pXiX2X3] = 
{X1X2X3) for the 3-measurement case. Hence, by using 
the different number of controlled gates in appropriate 
order we can calculate all the moments. The experimen- 
tally obtained moments are shown in Fig. 4. 

These experimentally obtained moments are inverted 
according to Eq. (4) to calculate the TTJP and are plot- 
ted along with the directly obtained TTJP using circuit 
shown in Fig. 1 as symbols in Fig. 5. The theoretical 
values for TTJP from moments and the one directly ob- 
tained are plotted as solid and dashed lines respectively. 
The results agree with the predictions of Eqs. (18) and 
(19) that the TTJP obtained directly and the one ob- 



tained from the inversion of moments do not agree. 



V. CONCLUSION 

In classical probability setting, statistical moments as- 
sociated with dichotomic random variables determine the 
probabilities uniquely. When the same issue is explored 
in the quantum context - with random variables re- 
placed by Hermitian observables (which are in general 
non-commuting) and the statistical outcomes of observ- 
ables in sequential measurements are considered - it is 
shown that the joint probabilities do not agree with the 
ones inverted from the moments. This is explicitly il- 
lustrated by considering sequential measurements of a 
dynamical variable at three different times in the specific 
example of a spin-1/2 system. An experimental investi- 
gation based on NMR methods, where moments and the 



joint probabilities are extracted independently, demon- 
strates the moment indeterminacy of probabilities, con- 
cordant with theoretical observations. 

The failure to revert joint probability distribution from 
its moments points towards its inherent incompatibility 
with the family of all marginals. In turn, the moment 
indeterminacy reveals the absence of a legitimate joint 
probability distribution compatible with the set of all 
marginal distributions - a common underpinning of var- 
ious no-go theorems in the foundational aspects of quan- 
tum theory. 
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